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Abstract. The BL Lac object Mkn 501 went into a very high 
state of activity during 1997, both in VHE 7-rays and X-rays. 
We present here results from observations at energies above 
250 GcV carried out between March and October 1997 with 
the CAT Cerenkov imaging Telescope. 

The average differential spectrum between 330 GeV and 
13 TeV shows significant curvature and is well represented by 
cf> E-£ +f3los ™ E ™\wiih: 

0o = 5.19 ± 0.13 stat ± 0.12 s y s - MC 

+1.66^-^ xl0 _H cm -2 s -l TeV -l ! 

a = 2.24 ± 0.04 stat ± 0.05 sys , and 

[3 = 0.50 ± 0.07 stat (negligible systematics). 

The TeV spectral energy distribution of Mkn 501 clearly 
peaks in the range 500 GeV-1 TeV. Investigation of spec- 
tral variations shows a significant hardness-intensity correla- 
tion with no measurable effect on the curvature. This can be 
described as an increase of the peak TeV emission energy with 
intensity. 

Simultaneous and quasi-simultaneous CAT VHE 7-ray and 
BeppoS AX hard X-ray detections for the highest recorded flare 
on 16 th April and for lower-activity states of the same period 
show correlated variability with a higher luminosity in X-rays 
than in 7-rays. 
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The observed spectral energy distribution and the corre- 
lated variability between X-rays and 7-rays, both in amplitude 
and in hardening of spectra, favour a two-component emis- 
sion scheme where the low and high energy components are 
attributed to synchrotron and inverse Compton (IC) radiation, 
respectively. 

Key words: : Galaxies: active - Galaxies: nuclei - BL Lacertae 
objects: individual: Mkn 501 - Gamma-rays: observations 

This article is dedicated to the memory of Eric Pare, one of the 
originators of CAT, who died in a road accident in July, 1998. 



1. Introduction 

Mkn 501, at z = 0.034, is the second-closest BL Lac object 
and, with Mkn 421 (z = 0.031), is one of the two extra-galactic 
sources thus far detected and confirmed at TeV energies. Its 
TeV emission was discovered in 1995 by the Whipple obser- 
vatory group ( Quinnet al. 1996 ) above 300 GeV, at the level 
of 8% of the TeV flux of the Crab Neb ula, and subsequentl y 
confirmed by the HEGRA collaboration ( Bradbury et al. 1997 ). 
Initially, Mkn 501 had been placed in the radio-selected BL 

from 1 Jy BL Lac and S4 ra- 



Lac family (Stickel et al. 1991 



dio samples). However, its spectral energy distribution (SED) 
in recent years resembled typical High-Frequency BL Lac ob- 
jects (mostly X-ray selected) with the low-energy peak in the 
EUV/soft X-ray energy band, making it a bright Einstein Slew 



Survey source (Ciliegi et al. 1995). In fact, it can be consid 



ered as being inter mediate between the rad io- and X-ray bright 
classes of BL Lac ( Lamer & Wagner 1998 ). 
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Fig. 1. Mkn 501 nightly integral flux levels above 250 GeV, expressed in units of 10~ n cm _2 s _1 , for observations in 1997. The 
average spectral shape as derived in Sect. 4 has been assumed in order to estimate the integral flux for observations far from the 
Zenith. 



From March to October 1997, Mkn 501 showed dramatic 
activity in the VHE 7-ray band, exhibiting flares during which 
it becam e up to ~ 8 times b righter than the Crab Nebula 
(e.g., see Protheroe et al. 1997); i.e., two orders of magnitude 



brighter than its discovery level. 

During the same period, Mkn 501 was also extremely active 
in the X-ray band as reported by the All Sky Monitor (ASM) on 
board the Rossi X-Ray Timing Explorer (RXTE). In April, high- 
intensity flares were recorded in the hard X-ray range by RXTE 
( |Lamer & Wagner 1998|), B eppoSAX ( JPian et al. 1998J ), and 
OSSE ( Catanese et al. 1997) on board the Compton Gamma 
Ray Observatory (CGRO). The synchrotron peak underwent a 
shift towards high energies as compared to archival data, shift- 
ing to energies of 100 keV. This unprecedentedly high energy 
made Mkn 501 the most extreme known BL Lac object in the 
X-ray sky. 

This paper is dedicated to the study of detailed spectral 
properties of this source, based on observations carried out with 
the CAT Cerenkov imaging telescope above 250 GeV in 1997. 

The next section describes briefly the CAT telescope and 
its specific event reconstruction method. The data-sample used 
for this paper and the resulting source light-curve are pre- 
sented in Sect. 3. Spectra derived for different-activity states 
are given in Sect. 4 and intensity-hardness correlations are stud- 
ied in Sect. 5. In Sect. 6 we compare simultaneous and quasi- 
simultaneous CAT and BeppoSAX data for the highest flare 
recorded on April 16 th to the lower activity states of the 7 th 
and 11 th of the same period. In Sect. 7 we discuss results and 
their implications for proposed emission models. The conclu- 
sions are given in Sect. 8. 



2. CAT Cerenkov Imaging Telescope 

The 17.8 m 2 CAT (Cerenkov Array at Themis) imaging tele- 
scope started operation in Autumn 1996 on the site of the for- 
mer solar plant Themis (Franc e). A complete descri ption of the 
CAT telescope can be found in Barrauet al. (1998) . 



The instrument records Cerenkov light from the particles in 
cosmic-ray showers. The analysis of th e shape and the pointing 
of the resulting images, as discussed in Le Bohec et al. (1998)[ 
takes full advantage of the very-high-definition camera. This 
latter has a 4.8° full field of view which is comprised of a cen- 
tral region of 546 phototubes in a hexagonal matrix spaced by 
0.13° and 54 surrounding tubes in two "guard rings". 

The 7-ray image analysis is based on the comparison of in- 
dividual events and theoretical average images as a function of 
impact parameter and energy. Accurate analysis of the longitu- 
dinal light profile and lateral width of the shower image is pos- 
sible thanks to the high-resolution camera. This fitting proce- 
dure allows the direction and energy of each 7-ray to be deter- 
mined with good accuracy: the angular resolution (0.14°) is of 
the order of the pixel size; when selecting showers with an im- 
pact parameter less than 130 m, the energy resolution is about 
20%, independent of energy, from approximately 300 GeV to 
15 TeV. At the highest energies, the image tails are cut by the 
field of view, but the light-profile analysis still permits a good 
energy estimation. 

It should be noted that, in general, for the Atmospheric 
Cerenkov Technique, there is a systematic uncertainty due to 
variations in the transparency of the atmosphere. Additional 
uncertainties due to unmastered detector drifts exist, but are 
fairly low for CAT, and are monitored within the data (e.g., us- 
ing muon ring images). A conservative estimate of the overall 
uncertainty on the absolute energy scale is 20%. 

To separate 7-ray initiated shower images from those due 
to hadronic showers we use the x 2 -like goodness-of-fit param- 
eter given by the above procedure, requiring the probability 
P x 2 > 0.35, and the pointing angle (which is the angle at 
the image barycentre between the actual source position and 
the reconstructed angular origin of the 7 ray), a < 6°. In addi- 
tion, only images with more than 30 photo-electrons and whose 
fourth brightest pixel has more than 3 photo-electrons are se- 
lected, in order to eliminate the lowest-energy triggers. These 
cuts yield a rejection factor of about 200, with a 7-ray effi- 
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ciency better than 40%. Further details on this specific method A = — 2 x log w^°\ - It behaves (asymptotically) like a x 



can be found in Le Bohec et al. (1998) 



distribution with one degree of freedom. 



3. Data Sample 

A total of 95.3 hours of data were taken on Mkn 501 (ON), to- 
gether with 36 hours on control regions (OFF), between March 
and October 1997. The data used for the analysis were selected 
by quantitative criteria including clear weather and stable de- 
tector operation. The data-set has been further limited to zenith 
angles < 44°, for which the detector calibration studies have 
been completed. This leaves 57.2 h ON and 22. 5 h OFF data for 
the analysis presented in the following sections. Fig. [I] shows 
the nightly flux levels as sampled by CAT above 250 GeV be- 
tween March and October, 1997. The average flux was approx- 
imately twice that of the Crab Nebula, while the highest flare, 
recorded on April 16 th (7-ray rate ^16 7/min), reached about 
8 times the Crab level with signal-to-background ratios of 2.7 
for the flare and 0.68 on average (see Fig. g(a,b)). 

4. Spectra 

4.1. Analysis Method 

As stated in Sect. 3, observations consist of ON and OFF source 
runs, the latter being used to estimate the background due to 
hadronic showers. Since the threshold energy of the telescope 
increases with zenith angle, ON and OFF data are classified 
within zenith angle bins (with a width of 0.05 in cos Z ). En- 
ergy bins have been defined so as to be larger than the tele- 
scope's energy resolution. Trigger and cut efficiencies - yield- 
ing the effective detection area A c s(E) - as well as the energy- 
resolution function are determined by detailed Monte-Carlo 
simulations of the telescope response as a function of energy 
and zenith angle. These simulations have been checked and cal- 
ibrated on the basis of several observables using the nearly-pure 
7-ray signal from the highest flare and muon rings. This will be 
reported in a forthcoming paper. 

The spectra are derived in two steps. First, the number of 
events passing the cuts within each energy and zenith-angle 
band is determined for ON and OFF runs. Then, for a given 
hypothesis on the spectral shape, a maximum-likelihood esti- 
mation of the spectral parameters is performed, taking into ac- 
count the effective area and the energy-resolution function of 
the telescope. 

Two hypotheses have been assumed for the spectral shape: 
a simple power law (Hq) and a curved shape (Hi). In fact, as 
we will see below, the data are not compatible with the for- 
mer shape. The curved shape is suggested by general consider- 
ations on emission processes within blazar jets (see Sect. 7.1). 
Whipple Observatory has also reported a curved shape in 
terms of a quadratic law in a log <j> vs. log E representation 
(Samuelson et al. 1998). We will use the same framework in 



4.2. Results 

In order to investigate spectral variations as a function of source 
intensity, three independent data-subsets corresponding to dif- 
ferent activity states of Mkn 501 - high-intensity flares (HF, 
integral flux > 53 x 10~ n cm _2 s _1 above 250 GeV), low- 
intensity runs (LF, integral flux < 12 x 10~ n cm~ 2 s _1 ), and 
mid-intensity flares (MF) - have been established. The aver- 
age over the complete data-set (AV) has also been considered. 
The quadratic law used to estimate the spectral parameters is 



b Q E T 



-(a+0log lo E To v) 



Results for each data-subset are given in 
Table [lj. Only statistical errors are quoted in this table. 

As noted in Sect. 2, the systematic errors on the absolute 
energy scale, conservatively estimated at 20%, lead to an error 
on the absolute flux depending on the spectral index (referred to 
below as "sys-atm"). Additional errors coming from the limited 
statistics of the simulated 7-rays (referred to below as "sys- 
MC") used to determine A c a{E) are twofold: 

- an overall scaling error which affects the absolute flux. 

- possible distortions of A g(E) as a function of energy 
within the Monte Carlo errors affecting the parameter a. 
The curvature parameter (3 is, however, almost unaffected. 

Both parameters have proved to be insensitive to modifications 
of the analysis procedure (which modify efficiencies and en- 
ergy resolution functions), showing that they are robust. 

For the AV spectrum, the complete results are: 

0o = 5.19 ± 0.13 stat ± 0.12 sys - MC 

+ 1 . 66 sys-atm 
-1.04 XlU 



-11 



L TeV~ 



order to allow for simple comparisons. 

To give an estimate of the relevance of Ho with respect to 
Hi, the likelihood ratio of the two hypotheses is defined as 



a = 2.24 ± 0.04 stat ± 0.05 sys , and 

= 0.50 ± 0.07 stat (negligible systematics). 

The corresponding covariance matrix for the fit, based on 
the statistical errors, is: 

0.01767 -0.00052 0.00443 \ 
-0.00052 0.00168 -0.00203 , 
0.00443 -0.00203 0.00475/ 
with a and f3 being correlated. 

The likelihood ratio, A, for Ho and Hi shows that the HF 
and MF states (and the AV spectrum) have significantly curved 
shapes: A has values of 29.1 and 47.1 respectively, correspond- 
ing to chance probabilities of 7 x 10~ 8 and 7 x 10~ 12 of ac- 
cepting Hi were the power-law hypothesis Ho true. For the LF 
state, the curvature significance is somewhat lower (a probabil- 
ity of 10~ 3 ), but this could be attributed to the weak intensity 
of the source. 

The SEDs are shown in Fig. [jf as can be seen, the VHE 
peak emission of Mkn 501 clearly takes place in the range 
of several hundred GeV. The corresponding maximum en- 
ergy, -E max , of the equivalent quadratic form - proportional to 
^Tev(^Tcv/S max )-' 31o g 1 o(^ v / £ — ) - is given by E max = 
10^ (see Table 0). 

There is some indication of an increase in E mlix with source 
intensity, which is equivalently shown by the variation of the 
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Fig. 2. The Mkn 501 signal seen with 
the CAT imaging telescope: a), b) a 
plots for April 16 th and for all data 
respectively (ON data shown by solid 
lines, OFF by points with error bars, 
ON-OFF inset). The signal is clearly 
seen in the direction of the source 
(small a), though the direction of 
some 7-rays is mis-identified, giving 
a small signal at a ~ 180°. c), d) two 
representations of the distributions of 
the reconstructed angular origins for 
all data (the scales in both plots are 
identical, the bin size is 0.05°, no 
background subtraction has been per- 
formed). The cross in c) marks the po- 
sition of Mkn 501. Events are selected 
by the P x 2 (i.e., shape) cut only. 



^TcV — 



cor- 



Table 1. Data-subset definition as a function of source intensity. The differential flux is given as d$/di?i 
f> £^ + ^ logl °- ETev) in units of 10~ n cm- 2 s- 1 TeV -1 . The fit has been done from 330 GeV to 13 TcV. $f and 
respond to a simple power-law hypothesis (Ho). A is the likelihood ratio of the two hypotheses Ho and Hi . The last two columns 

= 10^ 



give the peak-emission energies E max = 
obtained by re-fitting with j3 fixed at 0.5. 



(in TeV) for the second hypothesis, where for the final column, the values were 



Set 


ON 

(h) 


<pf 


a pl 


00 


Q 





A 


-C'max 

(GeV) 


p/3=0.5 

(GeV) 


LF 

MF 
HF 


13.6 

40.5 

3.1 


2.72 ±0.13 

4.10 ±0.10 

14.5 ±0.50 


2.45 ± 0.05 

2.46 ±0.03 
2.21 ±0.03 


3.13 ±0.19 

4.72 ±0.14 

17.6 ±0.61 


2.32 ±0.09 
2.25 ±0.05 
2.07 ±0.04 


0.41 ±0.17 
0.52 ±0.08 
0.45 ± 0.09 


10.7 
47.1 
29.1 


410 ±201 
583 ± 104 
840 ± 108 


521 ± 82 
563 ± 42 
890 ± 77 


AV 


57.2 


4.56 ±0.10 


2.46 ± 0.02 


5.19 ±0.13 


2.24 ± 0.04 


0.50 ±0.07 


61.5 


578 ± 98 


589 ± 38 



fitted values of a. Although the power-law hypothesis is clearly 
not favoured here, the corresponding variation of a pl values 
follows the same tendency. On the other hand, the curvature 
(given by (3) is compatible with being constant for the three 
data-subsets. Fixing j3 = 0.50 yields E max with smaller errors, 
as shown in the last column of Table [y. This reinforces a pos- 
sible correlation between intensity and E max . We will investi- 
gate this question in more detail in the next section by looking 
at intensity-hardness correlations. 



In order to determine the highest-energy 7-rays measured, 
the effects of the energy-resolution function and the uncertainty 
in the absolute energy scale must be taken into account. There 
is a signal at the 2cr significance level above 12 TeV, after cor- 
rection for the number of lower-energy events which migrate 
above this energy due to the energy resolution (and also to the 
lesser number which migrate to lower energy). Taking into ac- 
count also the uncertainty in the absolute energy scale implies 
(conservatively) that there is a 2a signal above 10 TeV. 
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Table 2. Hardness-ratio values for five different flux bands (zenith angle < 25°). The average integral flux per band above 
250 GeV is given in units of 1CP 11 cm _2 s _1 and has been computed assuming the average spectral shape. Set I corresponds to 
the highest flare observed on April 16 th and set II to the second-highest flares of April 13 th and June 7 th . 



Set ON (h) 



Flux 



R 



[>900/>450] 



R\ 



>1200/>600] 



Rt 



>1500/>900] 



1 


1.5 


68-105 


0.56 ± 0.02 


0.50 ±0.02 


0.56 ± 0.03 


II 


3.5 


45-68 


0.53 ± 0.02 


0.49 ±0.02 


0.56 ± 0.02 


111 


5.5 


28-45 


0.47 ± 0.02 


0.44 ± 0.02 


0.53 ±0.03 


IV 


10.5 


17-28 


0.48 ± 0.02 


0.42 ± 0.02 


0.49 ± 0.02 


V 


20.5 


5.5-17 


0.45 ± 0.02 


0.41 ±0.02 


0.51 ±0.03 



450-900 GeV 



600-1200 GeV 



1 ' 1 ' 1 


900- 1500 GeV 



Fig. 4. Hardness-ratio values 
vs. source intensity for three 
different energy bands. Inten- 
sities are given as the integral 
flux above 250 GeV in units 
of 10 _u cm^s -1 . 



5. Spectral Hardness vs. Source Intensity 

As compared to a usual spectrum analysis method, a more ro- 
bust manner to check for spectral changes is to consider two 
given energies E\ ov/ and E m id, and the hardness ratio defined 



as: R = 



iV E ; 



JVe; 



Given the large energy-bins, this ratio is rather insensitive 
to statistical fluctuations and entails the estimation of only one 
parameter. Nevertheless, as the detection threshold of atmo- 
spheric Cerenkov telescopes increases with zenith angle and 
is sensitive to slight changes in weather conditions, it is desir- 
able to limit the data-set to low zenith angle runs (here 25°) 
and choose E\ ov/ well above the telescope trigger threshold of 
250 GeV (here 450 GeV or above). 

Although these two restrictions reduce the available statis- 
tics, they allow for minor acceptance corrections by weighting 
data to normalize them to those taken at the Zenith. This mini- 
mizes any possible systematic effects. Furthermore, to achieve 
the best energy resolution, an additional cut on shower impact 
parameter (< 130 m at the Zenith) has been introduced. 

The limited statistics still forbid a reliable computation of 
the hardness ratio on a run-by-run (or even night-by-night) 
basis, except for the highest-intensity runs. Therefore, data 
have been divided into five sets with different average fluxes. 
For each data-set in flux, the hardness ratio has been com- 
puted for three different energy bands, [E\ ow , -Emid]GeV = 
{[450, 900]; [600, 1200]; [900, 1500]}, hereafter referred to as 

#[>900/>450] . -R[>i200/>600]' and -R[>i500/>900] ■ The number 
of hours for each data-set and the corresponding average flux 
and hardness ratio values are given in Table g. 



As can be seen in Fig. 0, a clear correlation between in- 
tensity and hardness is present for all energy-band definitions. 
The chance probabilities of observing such an effect in the case 
of an absence of correlation for i?[> 90 o/>450], #[>i200/>600] 
and i?[>i5oo/>90o] are respectively: 3.8 x 10~ 4 , 9 x 10~ 4 , and 
7 x 10 . We note that as E\ ow and E mi( \ are raised to higher 
energies, the significance of the correlation decreases owing to 
lack of statistics. 

Considering the curved spectral shape measured in the pre- 
vious section, and fixing the curvature at (3 — 0.50, we can cal- 
culate the shift in E max necessary to cause the observed max- 
imum change in R. The variation is of the order of 400 GeV, 
which is fully compatible with the results of the previous sec- 
tion seen in the final column of Table |l| 

6. Hard X-ray and 7-ray correlation 

During the April 1997 multi- wavelength campaign, Mkn 501 
was observed from X-ray to VHE energies (e.g., see 
|Catanese et al. 1997 ). In the energy range 0.1-200 keV, the 
BeppoSAX satellite observed for three periods of ^10 hours: 
on April 7 th , 11 th and 16 th (|Pian et al. 1998|). The RXTE satel- 



lite observed in the 2-80 keV range ( Lamer & Wagner 1998J ) 
between July 11 th and 16 th . X-ray and high energy 7-ray ob- 
servations were also performed by the OSSE and EGRET in- 
struments respectively, on board CGRO. In the VHE regime, 
CAT, HEGRA, Whipple and other instruments were active 



(Protheroeetal. 1997) 



In the X-ray range, extremely hard spectra were reported 
by BeppoSAX indicating that the synchrotron power output 
peaked at 10 keV or more. OSS E observations gave similar 
conclusions (Catanese et al. 1997). 
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Table 3. Observation logs and power output (in units of 10 -11 erg • cm~ 2 s _1 ) for BeppoSAX ( ^ian et al. 1998 ) and CAT. For 
CAT, the integral fluxes corresponding to April 7 th and 11 th are calculated assuming the spectral shape from the LF and MF 
data-subsets, respectively. 



Date 



BeppoSAX 
Observing time (UT) 



X-ray power output 
13-200 keV 



CAT 

Observing time (UT) 



VHE power output 
0.25-15 TeV 



April 7 th 
April 11* 
April 16* 



05:30: 
05:53: 
03:36: 



29^ 16:01:47 
57 -> 16:25:49 
52 -> 14:35:48 



37.5 ± 1.5 

51.5 ±1.5 

158.0 ±2.0 



00:34:06 -> 03:56:53 
01:06:57 -+ 03:53:16 
01:42:34 ->■ 03:45:55 



9.8 ±2.4 

22.4 ±3.1 

100.0 ± 7.5 



E(TeV) 



1th 



-\th 



High— intensity Flare (HF) 



Mid— intensity Flare (MF) 





t 

-10 J» 



■'(Hz) 

Fig. 3. vF v (or E 2 d$/dE) spectra for the HF, MF, and LF 
states. Curves represent the global flux estimated by method 
described in Sect. 4.1. The points, corresponding to the indi- 
vidual bin intensities, are only indicative and were obtained 
by the following procedure: The expected number, JVfi t , of 7- 
ray events in each energy bin is computed using the fitted flux 
parameters and taking into account the acceptance and energy 
resolution functions of the telescope. The observed number of 
7-ray events, Nobs, is obtained by subtracting ON and OFF data. 
The ratios -j^ give a good approximation of the fluctuation in 
each bin, and are used to correct the flux values yielded by the 
fitted curve. Error bars give the combination of uncertainties on 
iVobs and on the effective detection area (statistical only). 



The highest-intensity TeV flare recorded by CAT in 1997 
occurred on April 16 th (see Fig. nh, where CAT and BeppoSAX 
observations overlapped for about ten minutes, providing some 
simultaneous VHE 7-ray and X-ray data. The observation logs 
are given in Table 0. On April 7 th and 11 th , the source was 
less bright and only quasi-simultaneous data were available, 
as there was a delay of about two hours between CAT and 
BeppoSAX pointing periods. However, data from all Cerenkov 
telescopes (CAT, Whipple and HEGRA) show no large varia- 



tions in source flux for the periods 5 to 9 and 10 to 12 



)th 



(Protheroeet al. 1997). Hence, 7-ray and X-ray levels can be 



considered as simultaneous data and compared in a reliable 
manner. 

Table g shows comparisons of the VHE and X-ray power 
output from BeppoSAX and CAT data. The ratio of VHE to X- 
ray power outputs is plotted in Fig. gas a function of the X-ray 
power, and clearly increases with source brightness. 

The spectra are shown in Fig. g. Because of the very high 
intensity of the source on April 16 th , it is possible to derive 
its spectrum for that single observation (one and a half hours of 
data). In order to represent the spectrum corresponding to April 
7 th , the LF data-subset comprising runs with similar flux levels 
(within 30%) was used. 

Four characteristics come out of these comparisons: 

(i) the fluxes and SED shapes seen in the CAT and BeppoSAX 
data, together with the EGRET upper limit in the GeV 
domain, imply a double structure in the power output of 
Mkn 501; 

(ii) hard X-rays and VHE 7-rays show correlated variability; 

(iii) the luminosity in X-rays remains higher than in 7-rays; 

(iv) during the flares, the VHE intensity recorded by CAT in- 
creased by a larger factor than the X-ray intensity recorded 
by BeppoSAX (Table |). Therefore, the ratio of VHE to X- 
ray intensities increased with source brightness (Fig. ph. 




100 150 

X— ray power output (1 0"' 



Fig. 5. The ratio of VHE to X-ray power outputs as a func- 
tion of the X-ray output, from the data in Table g. On this plot, 
a linear relation between the VHE and X-ray intensity would 
appear as a horizontal line, and a quadratic relation as a line 
passing through the origin. 



A. Djannati-Ataiet al.: VHE Spectral properties of Mkn 501 with CAT in 1997 



10 10 10 10 10 



E (MeV) 

10 7 



ieppoi 



0&\ 

Li\4 



16 April 



-10 .«• 



***u 



1 



r *m 



7 April 



EGRET 

T 



,(LF) 



f - 10 

\ : 



f 



\ 



+ .10 

'..I 

! : 

\ : 



Mrk 501 April 97 



(/(Hz) 

Fig. 6. X-ray and VHE spectra given as vF v . For April 7 th and 
16 th , BeppoSAX data are from Pian et al. (1998). The LF sub- 
set in CAT was used to represent the spectrum corresponding to 
April 7 th . The EGRET upper limit is taken from Samuelson et 
al. (1998) and corresponds to observations between April 9 th 
and 15 th . 



7. Discussion 

CAT observations of Mkn 501 during 1997 have revealed sev- 
eral spectral properties: a curved shape and an emission ex- 
tending above lOTeV (Sect. 4), hardness-intensity correlations 
(Sect. 5), X-ray and 7-ray correlations (Sect. 6). Together with 
the power deficit in the GeV range, this strongly suggests a 
two-component emission spectrum for Mkn 501. 

It is generally accepted that the high luminosity of blazars 
(~10 45-48 erg/s) requires a relativistically-amplified beamed- 
emission. Beaming is also required for the source to be trans- 
parent to 7-rays if both VHE and IR to UV emissions orig- 
inate in the same region of the jet (Dermer & Gehrels 1995; 
Dondi & Ghisellini 1995). Two major processes have been 
suggested as the dominant source of the high-energy emis- 
sion from blazars: inverse Compton (IC) scattering of low- 
energy photons by ultra-relativistic electrons, and 7r°-decay 
through photo-production or secondary production of pi- 
ons by an ultra-high-energy hadronic component in the jet 



( Mannheim & Biermann 1992; Mannheim 1993) 



In the IC case, a double-structure spectrum in terms 
of synchrotron and IC emissions is predicted with the two 
components being correlated in intensity, hardness, energy 
of the peak luminosities, etc. The correlations depend on 
the nature of the target-photon population. This may con- 
sist of either synchrotron photons within the jet - the syn- 



chrotron self-Compton , or SSC, mechanism (Konigl 1981; 
Marscher & Gear 1985 ghisellini et al. 1996) - or external 
radiation, either from the accretion disk ( |Sikora et al. 1994 ; 



Dermer et al. 1992), or reflected by the broad-line region 



clouds (Ghisellini & Madau 1996). Spectral variations of the 



primary electrons (due for example to a change in the injection- 
acceleration efficiency or in cooling conditions), seen through 
their effect on the synchrotron radiation, should be reflected in 
the high-energy emission, especially when the IC interactions 
take place in the Thomson limit. In the Klein-Nishina regime, 
hardening effects can be attenuated because of a severe cross- 
section drop at high energies. 

In hadronic models, no particular predictions are made for 
a marked two-component emission. A continuum emission is 
predicted due to secondary cascades produced by hadronic in- 
teractions of ultra-high-energy protons, which tend to fill the 
gap between the high-energy (7-ray) and low-energy (X-ray) 
domains. 

7.1. Curved Shape 

The spectral shape derived in Sect. 4 between 330 GeV 
and 13 TeV is clearly curved. This curvature is in contrast 
with reports on the other TeV blazar, Mkn 421, whose spec- 
t rum has been observed to be compati ble with a power-law 
( Krennrich et al. 1999 ; Piron et al. 1999 ). The curved shape is 
predicted by IC models of blazar spectra but can only be seen 
if the IC component is observed over a broad dynamic range 
close to the maximum power emission energy (i.e., close to the 
maximum of the IC bump in a vF v representation). The cur- 
vature could arise from intrinsic absorption effects (within the 
jet or due to nearby radiation fields, from the accretion disk, 
for example) and/or from a cutoff in the primary particle spec- 
trum. However, when observing within a limited energy range 
- as compared to the energy scale over which the emission 
phenomenon is taking place - not including the peak emission 
region, a power-law spectrum is naturally seen. We note that 
the reported synchrotron peak of Mkn 501 during 1997 is well 
above that of Mkn 421 (~100 keV as compared to ~1 keV) 
and that the IC and synchrotron peak-emission energies are 
expected to be correlated in IC schemes. This correlation has 
also been established on a purely observational basis using a 



large sample of blazars (Fossati et al. 1998). Thus, the differ- 



ence in observed spectral shapes for Mkn 421 and Mkn 501 
may well arise because the peak of the high-energy component 
for Mkn 421 is below the VHE energy domain, whereas the 
peak for Mkn 501 is well above our threshold. 

The overall spectrum, and particularly the highest energy 
7-rays (10-20 TeV), may be affected by intergalactic absorp- 
tion due to pair creation on the intergalactic background IR ra- 



diation ( Stecker&de Jager 1998 ; Billeretal. 1998 ). Therefore 



any interpretation of VHE spectra in terms of emission models 
should take this effect into account. It could be disentangled 
from intrinsic absorption features if we dispose of a large sam- 
ple of VHE emitters. Nevertheless, the different spectral shapes 
of Mkn 501 and Mkn 421 - which have similar red-shifts and 



have been detected above 1 TeV ( Krennrich et al. 1 999 ) - im 



ply already that the curvature of the former must be mainly due 
to intrinsic features. 
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7.2. Hardness-intensity and ^-ray/X-ray Correlations 

In Sect. 5, a correlation between VHE intensity and spectral 
hardness was demonstrated. If we assume a constant curvature 
for different-activity states as indicated by the data, this effect 
is equivalent to an increase of approximately 400 GeV for the 
peak emission energy E max for an intensity change of about an 
order of magnitude between the LF and HF subsets. 

BeppoSAX results have established that the energy of the 
maximum power output during the a ctive period is mu ch larger 
than during the low-emission state ( Pian et al. 1998 ). A com- 
parison of April 7 th and 16 th X-ray data suggests that this hard- 
ening continues as the intensity increases during the flare. 

This simultaneous observation of flaring activity and spec- 
trum hardening, in both X-ray and VHE energy ranges, is nat- 
urally explained by a common primary electron population, as 
expected in the framework of inverse Compton models. Pian et 
al. (1998) claim that a significant change in the magnetic field 
intensity or the Doppler factor is unlikely since the expected 
changes in the lower-energy part of Mkn 501 spectrum which 
would result have not been observed; i.e., the X-ray emission 
below 0.5keV did not show any dramatic variation. In this case, 
this would leave the injection of a transient electron population 
and/or a change in the acceleration efficiency (an increase of 
the maximum energy of electrons) as the most likely source of 
the flares. 

In the framework of IC models, the observed increase of 
the VHE to X-ray intensity ratio with source brightness (see 
Fig. ph is qualitatively in agreement with the SSC mechanism 
which predicts a quadratic increase of VHE intensity with X- 
ray brightness due to the fact that both the electron popula- 
tion and the "target" synchrotron populations are enhanced. 
The yield of the SSC contribution to the IC process should be 
damped by the decreasing cross-section in the extreme rela- 
tivistic Klein-Nishina regime. The highest intensity data point 
in Fig. q] is consistent with such an effect. 

7.3. Comparison with other VHE detections 

The Whipple group's results for Mkn 501 between 300 GeV 
and lOTeV (Samuelson et al. 1998) show an average spectrum 
with a = 2.20 ± 0.04 stat ± 0.05 syst and = 0.47 ± 0.07 stat , 
in good agreement with CAT results within the statistical er- 
rors. HEGRA ( Aliaronian et al. 1999) has also demonstrated a 



curved spectral shape. Given the source variability and the dif- 
ferent source sampling due to local conditions, one does not ex- 
pect the absolute flux levels observed in the three experiments 
to be identical. Nonetheless, the convergence of these three ob- 
servations as regards the spectral shape is noteworthy. 

No hardening effect has been claimed by either Whipple 
or HEGRA. While this appears to contradict the effect seen 
in our data, it may be explained for each case. We note that 
Whipple's reported results are limited to the period between 
February and June 1997 and include 21.1" of data. The CAT 
results are based on a larger sample (40. 5 h of observations for 
the hardness-intensity studies), and include runs from June to 



October 1997, providing large statistics of the source emission 
in a lower-activity state. The effect reported here could hardly 
be seen with lower statistics. 

Considering HEGRA results, their higher energy threshold 
might well result in a lower sensitivity to the hardness-intensity 
correlation: as shown in Sect. 6, this effect becomes less signif- 
icant when we consider higher energy events. 

8. Conclusion 

The results reported in this paper are based on the large sample 
of data taken by CAT from March to October 1997, includ- 
ing periods with very different source intensities. The Mkn 501 
spectrum has been measured from 330 GeV to 13 TeV. The 
curvature reported by the Whipple group is confirmed. In ad- 
dition, CAT data show an intensity-hardness correlation which 
can be simply described by a shift of the peak TeV emission 
energy. 

The overall picture emerging from these spectral proper- 
ties and from the comparison with simultaneous and quasi- 
simultaneous BeppoSAX observations favours the interpreta- 
tion in terms of an energetic electron beam propagating in a 
magnetic plasma and producing X-rays through synchrotron 
radiation as well as VHE 7-rays through inverse Compton scat- 
tering of low-energy photons. 
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